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Abstract. We discuss expectations of saturation physics for various observables in heavy ion 
collisions. We show how simple saturation-inspired assumptions about particle production in heavy 
■ ion collisions lead to Kharzeev-Levin-Nardi model. Comparing this model to RHIC data on particle 

^jrj| multiplicities we conclude that saturation effects may play an important role in particle production 

and dynamics at the early stages of Au—Au collisions already at RHIC energies. We then estimate 
the contribution of the initial state two-particle azimuthal correlations to elliptic flow observable V2 
in Au —Au collisions by constructing a lower bound on these non-flow effects based on v? obtained 
^sO . from the analysis of proton-proton (pp) collisions. 

! INTRODUCTION 

r- ■ 
o ■ 

Saturation/Color Glass Condensate physics is based on the observation that the small-* 
wave functions of ultrarelativistic hadrons and nuclei are characterized by a hard scale 
Q s , known as the saturation scale JESS 01- The scale Q s arises due to saturation 
of partonic densities at small-x and is an increasing function of energy and atomic 
number of the nucleus q, This large scale makes the strong coupling constant 

small 0Cy(<2t) <C 1 leading to dominance of the classical gluonic fields in all high energy 
processes \4aM - Gluon production in high energy collisions is given by the classical field 
of the scattering color charges [8]. Corresponding gluon production cross section was 
found for pA collisions in |9|] and the effects of quantum evolution ||5|] were included in it 
in [10]. The gluon production cross section for heavy ion collisions (AA) at the classical 
level has been studied both numerically ifllll and analytically [12]. Since it is quite not 
clear at present how to include the effects of nonlinear quantum evolution o] in the 
results of lflll[l2ll . one has to construct models to describe the actual rapidity-dependent 
data produced in heavy ion collisions. Below we are going to show how some of these 
models, based on rather basic properties of saturation physics, provide a reasonably good 
description of RHIC data. 



PARTICLE MULTIPLICITY FROM SATURATION MODELS 

Multiplicity at Mid-Rapidity Versus Centrality 

Classical field ~ 1/g leads to produced gluon multiplicity 

dN /A A \ 1 
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Gluon transverse momentum spectrum described by a single scale Q s can be written as 



with some unknown function f(kj_/Q s ) to be determined by actual calculations. Inte- 
grating over k and b yields 

dN 1 . 2 

— = const — TlR 2 Q 2 (3) 
ay a s 

where the value of the constant is determined from f(k±/Q s ). Following lfl3l fbUl we 
assume that the scale for the coupling constant in Eq. © is set by Q s . (This step, of 
course, goes beyond the classical limit and assumes that Eq. © is valid even when 
running coupling corrections are included.) Then, as R 2 ~ A 2 / 3 and if Q 2 ~ A 1//3 0,0], 
together with running coupling 

we conclude from Eq. © that 

For heavy ion experiments at different collision centralities we substitute A by the 
number of participants N part so that Eq. © becomes 

1 dN 



Eq. © allowed the authors of 11311 to correctly predict the particle multiplicity at mid- 
rapidity at RHIC as a function of centrality at y/s = 130 GeV. A fit of particle multiplicity 
at other values of rapidity at y/s = 130 GeV taken from [14] is shown in Fig.[H 
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FIGURE 1. Saturation model fit of the PHOBOS data on total charged particle multiplicity at mid- 
rapidity as a function of centrality at y/s = 130 GeV taken from 1 14]. 



Multiplicity as a Function of Energy 



Eq. © allows one to test whether the scaling of total particle multiplicity with energy 
is consistent with saturation/Color Glass predictions. Using the fact that Q% ~ 1/xgj 
10,0] in Eq. © and, for the moment, dropping the slower (^-dependence in a s leads to 



dN / df]{y/si) Vv 7 ^ 
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Using PHOBOS data for total charge multiplicity at y/s =130 GeV for most central 
collisions 

dN 

— (y/s = 130GeV) = 555 ±\2( st at) ± 35 (syst) (8) 
dr\ 

together with A = 0.25 -=-0.3 obtained in lfl5ll by analyzing HERA data, Kharzeev and 
Levin ifbUl predicted the total charge multiplicity at y/s = 200 GeV to be 

dN 

— ( v / s = 200GeV) = 616^634 (9) 
dr\ 

in agreement with the later measured PHOBOS result 
dN 

— (v^ = 200GeV) = 650±35{stat). (10) 



^ Versus T] and N par t 

Describing rapidity distribution of the produced particles requires a little more mod- 
eling. Assuming ^-factorization for particle production cross section 




FIGURE 2. Saturation model fit of the PHOBOS data on total charged particle multiplicity as a function 
of rapidity and centrality at y/s = 130 GeV taken from 1 14]. 
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along with saturation-inspired unintegrated gluon distribution functions (fa(k) ~ a s/t? 
if > Q s and fa(k) ~ S±/ CC s if fc^ < 2^) the authors of ifbUl produced an impressive 
fit of the PHOBOS data on charged particle multiplicities as functions of rapidity and 
centrality at y/s =130 GeV shown in Fig.|2j The predictions made in Ill4ll for similar 
multiplicity data at y/s = 200 GeV were also in good agreement with the later published 
BRAHMS data. 

Phenomenological success of the saturation models presented above does not contra- 
dict the possibility of strong final state interactions leading to formation of quark-gluon 
plasma. As was argued in IU6I1 . thermalization in the saturation framework would not 
introduce any fundamentally new scales leaving Eqs. © and © practically unchanged. 
Late stage interactions are also not very likely to significantly modify the rapidity distri- 
bution of Fig.Eldue to causality constraints. 



NON-FLOW CONTRIBUTION TO V 2 

The contribution of non-flow two-particle azimuthal correlations from the early stages 
of heavy ion collisions to the elliptic flow observable V2 has been estimated in Ill7ll 
using a saturation-inspired model of particle correlations. Here we are going to construct 
a model-independent lower bound on these non-flow effects using V2 extracted from 
the analysis of pp data. We start with the definition of V2(pr) for pp corrected for 
uncertainty in the "reaction plane" definition (or, equivalently, defined through the two- 
particle correlation functions, such that V2(pr) <V2> = < cos2(0 p — fa) >k ) 
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If we assume that the relative magnitude of the correlated terms in Eq. (fT2l compared to 
uncorrelated ones is roughly the same for all px, we can drop the former compared to 
the latter finding that Eq. (fT2b is approximately bounded from above by 
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We want to estimate the contribution to v^ipr) of the non-flow correlations of the same 
physical origin as the ones giving rise to v^ipr) m Eq. (fT2l . The contribution is 
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Using the fact that, approximately, both in saturation models and in the data dN/dy 
Npart lll3Lll4ll . we rewrite the right hand side of Eq. ([T3l as 



/ AjAA 

vf(/>r) < vf WU-/ to ]J (15) 

with and N pp total particle multiplicities in AA and pp collisions, proportional to 
the average number of participants involved. Inverting Eq. (fl3t we obtain a lower bound 



/ NPP 

V^(PT)\non-flo W > vf (/>_ ) \J J^- (16) 

Preliminary analysis of pp data yields v^ipr) ~ 1 athigh-py [18]. To get a lower bound 
we use N^ rt = 394 and Nf art = 2 in Eq. (__]> obtaining vf (p T ) \ non - f i ow > 7%. 
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